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Abstract. We examine the propagation of collisionless particles emitted from a 
spherical shell to infinity. The number distribution at infinity, calculated as a function 
of the polar angle, exhibits a small deviation from uniformity. The number of particles 
moving from the polar region toward the equatorial plane is slightly larger than that 
of particles in the opposite direction, for an emission radius > 4.5Af in extreme Kerr 
space-time. This means that the black hole spin exerts an anti-collimation effect on the 
particles stream propagating along the rotation axis. We also confirm this property in 
the weak field limit. The quadrupole moment of the central object produces a force 
toward the equatorial plane. For a smaller emission radius r < 4.5M, the absorption of 
particles into the black hole, the non-uniformity and/or the anisotropy of the emission 
distribution become much more important. 



PACS numbers: 04.20.Cv, 04.70.-s, 95.30.Sf, 97.60.Lf, 98.38.Fs, 98.58.Fd 
1. Introduction 

Astrophysical relativistic jets are commonly observed in AGNs, microquasars and 
gamma ray bursts. The Lorentz factor of the outflow motion is 7 ~ 10 in AGNs, 
7 « 2 — 5 for microquasars and 7 ~ 100 for gamma ray bursts. It is considered likely that 
a central black hole is related to the formation of such jets due to the large luminosity 
and short timescale of the jets. However, the actual mechanism of such collimated 
outflows has not yet been resolved, although many theories have been proposed. One 
possible mechanism that has been suggested relates to magnetic hoop stress. Recently, 
large-scale MHD simulations of jets have been successfully performed by several authors. 
For example, McKinney[T] numerically studied the hydrodynamical flow emitted from 
an accretion disk, where the jet formation extends to a distance of ten thousand times 
the radius of the central black hole. In order for a collimated structure to form, the 
magnetic pressure has to be sufficiently large in comparison to the ambient pressure of 
the surrounding medium. 

In this paper, we explore an alternative possibility. That is, we investigate a 
mechanism involving the black hole spin itself. The strong gravitational field near a black 
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hole also affects the flow of particles(e.g.,[2J [3]). The effect of gravity may be important 
when the emission region is concentrated in the vicinity of the central black hole. Bicak 
et al. [I] discussed this effect by calculating the geodesies of some particles in Kerr space- 
time, de Felice and Carlotto [5] took into account the slow loss of angular momentum and 
energy of the particles moving along the geodesies. They derived a certain condition for 
such dissipative collimation and discussed its astrophysical relevance [6], [7]. Semerak[8] 
considered the effect of a spinning particle in Kerr space-time. More recently, Gariel 
et al.[9] discussed the collimation of geodesies by using cylindrical coordinates. These 
authors pointed out the possibility of collimated outflow occurring along the rotation 
axis of a Kerr black hole. The axis is a unique preferable direction of the flow, and 
some trajectories are in fact concentrated along the axis. Despite the fact that the 
possible collimation mechanisms in Kerr space-time have been discussed extensively, 
it is not clear what fraction of the flow (the ratio of the flow along the rotation axis 
to the overall emission) is collimated. For example, if particle emissions are selected 
with orbital parameters suitable for the flow along the axis, then the collimated flow 
is realized. It is therefore important to evaluate the degree of the collimation and the 
results will be useful for the astrophysical application. 

The problem is closely related to the emission process itself. To tackle this problem, 
we present a model for studying geometrical collimation. In our model, billions of 
particles are isotropically emitted from a shell with a constant radius covering a central 
black hole, and the trajectories are numerically calculated. We use a test-particle 
(dust) approximation to calculate each trajectory. The collision of particles is ignored 
in order to evaluate the geometrical effect in an idealized situation. The collisional 
hydrodynamics may work toward local isotropy, and may partially hide the effect. The 
angular distribution of the particles at the emission radius is assumed to be uniform, 
but at infinity some deviation from a uniform distribution is produced by Kerr black 
holes. In this way, we study to what extent a rotating black hole affects the collimation 
of collisionless particles. Our model is not rigorously realistic, but is instead an idealized 
mathematical model which is useful in evaluating what fraction of the flow is collimated. 

The organization of this paper is as follows. In Section [21 we summarize our models. 
The numerical distribution of the trajectories of the particles at infinity is calculated as a 
function of the polar angle, the particles having been emitted with a uniform distribution 
from a 'spherical' shell. It is not easy to define a spherically symmetric surface in 
non-spherical space-times. Two dimensional spatial surface with a constant radius is 
approximately regarded as a sphere. This definition, which is valid for asymptotically 
large radius, is used in this paper. The numerical results in Boyer-Lindquist coordinate 
are given in Section [3j In order to consider the physical mechanism, the same calculation 
is performed in weak gravity. Asymptotically Cartesian and mass centered to order N 
(ACMC-N) coordinate is used and the results are given in Section HI As the emission 
radius becomes smaller, a larger number of particles which are emitted outwardly are 
absorbed into the central black hole. In this case, the spatial uniformity and the isotropy 
of the particle distribution on the emission sphere become meaningless. In order to 
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demonstrate the absorption effect, some results are also shown in Section Finally, a 
discussion of our results is given in Section [61 We use units of c = G = 1 in this paper. 



2. Model and method 



The Kerr space-time for a black hole with mass M and angular momentum aM is 
described in the Boyer-Lindquist coordinates (t, r, 9, 0) as 



ds 4 
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There are three constants related to the geodesic motion in Kerr space-time: the energy 
E, the z component of the angular momentum L z and Carter's constant Q [TtH [TTj [T2]. 
The 4-momentum p^ of a particle with a rest mass m is therefore written as 



-E, ± A , 



ie 1 / 2 , l 5 



where 



i? = (E(r 2 + a 2 ) - L z a) 2 - A [m 2 r 2 + (L 2 - aE) 2 + Q] , 
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In order to allow for the emission model of particles to be considered, we introduce 
a reference frame (t, f , 9, <p) measured by zero-angular momentum observers (ZAMO) 
(e.g., [12]). The reference frame is represented by four orthogonal unit vectors, i.e., one 
timelike vector em) and three spacelike vectors e( n )(n = 1, 2, 3). The explicit components 
are given by = [e~ v , 0, 0, ue~ v ], ef 1} = [0, A 1/2 /p, 0, 0] , ef 2) = [0,0,l/p,0] and 



'(3) 



[0, 0, 0, e ^] . The 4-momentum p of a particle is expressed by 



p = -Ete(o) + p ( cos #e(!) + sin 9 cos 0e( 2 ) + sin 9 sin 0e( 3 ) 



(6) 



where the angles 9 and <fr are determined by the emission direction of the particle. The 
energy and the linear momentum satisfy the special relativistic relation p = (E 2 



2\l/2 



E = r?i7, where 7 
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and v is the velocity of the test particle. 



From Eqs.([3]) and ([6]), the constants of motion E, L z , Q can be expressed by 
quantities measured by a ZAMO, 
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We consider a particle moving along the geodesic curve without any other forces, 
so that r-6 motion can be calculated by integrating the differential equation 



with initial condition (r,9) = (r e ,8 c ). The numerical integration method is a fourth- 
order Runge-Kutta method with an adaptive step size, and the numerical errors are 
adjusted to be small enough. In order to check the validity of our numerical integration, 
we calculated the same problem as Bicak et al.[l] and obtained the same results. 

We consider the particle distribution, f(cos8), at infinity r^, when a large number 
of test particles, N tot , are emitted from a shell of radius r c . The particle distribution 
at r e is assumed to be uniform, that is, the number of particles dN emitted from 9 C to 
9 C + dO c is dN = (iV tot /47r)(27rsin 6 e d8 e ) for ^ 9 C ^ n. As for the emission model at 
(r e , 9 e ), we also assume that the particles have the same energy E, and that the emission 
direction is isotropic with respect to the angles 9 and <fi. The angle 9 might be limited 
to ^ 9 ^ 7r/2, since most of the particles satisfying the emission angle tt/2 ^ 9 ^ tt 
are absorbed by the black hole and can therefore be ignored from the beginning. In 
the numerical calculation, the directions 9 and <ft are randomly chosen at r c , and the 
constants of motion can be calculated by Eqs.(J7])-(j9]). The asymptotic value 9 of each 
particle is numerically evaluated by solving Eq.( jT0i) to r^, where r m is a large number 
satisfying 3> M, which is set to 'infinity' in the numerical calculation. By calculating 
iVtot trajectories, we obtain the number distribution f*(cos9). 

In order to study the geometrical effects, we are not interested in a fine-grained 
distribution /* but in a coarse-grained distribution, i.e., an averaged distribution on a 
larger angular scale. In the numerical calculations, we divide cos 0(0 ^ cos 9 ^ 1) into 
Nq segments with a small interval d cos 9, and count the number of the trajectories which 
arrive at the region between cos# and cos 9 + dcos9. The number is approximated by 
/ (cos 9)d cos 9, and the function f(cos9) is numerically constructed. 

In our model, the isotropy of the emission direction and the spatial uniformity are 
assigned by a ZAMO at a constant radius in Boyer-Lindquist coordinates. If the emission 
is uniform with respect to the proper surface element at r e , then the number of particles 
dN should be modified as dN = g (N tot /A.h){2ti sin 9 e d9 e ) for ^ 9 e ^ 7r. The correction 
factor g is estimated as g = 2A 1/2 / (J* A 1/2 sin 9d9) w l + a 2 (r -2M)(3 cos 2 9- l)/(6rf), 
which gives a small variation of 5 x 10~ 3 for r e = 4.5M. The correction \g — 1| seems to 
increase as r e decreases. We have numerically evaluated the factor and found that there 
is a maximum value at a certain radius r e . For example, in extreme Kerr space-time 
a = M, we found that \g — 1| < 1.5 x 10 -2 , where the maximum occurs at r e ~ 2.4M 
and 9 = 0, n. We actually include this correction factor into the numerical calculations 
in Section 3 and 5, and the obtained results are not significantly different, since other 
physical effects are much larger in the magnitude. 
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Figure 1. Deviation from uniform distribution as a function of cos 6. The left panel 
shows the results obtained for a flat space-time, and the right panel shows the results 
obtained for a Schwarzschild space-time. 



3. Effect of geometry on propagation 

In principle, the number distribution function f(cos8) ranges from to A^ to t and clearly 
depends on N tot . Therefore, this property is not suitable when comparing, for example, 
numerical results obtained using different total numbers of particles. The deviation 
from uniform distribution is more important. We therefore use the following normalized 
distribution F(cos6) defined by 

f(cosO) ^ /(cos g ~ W x 10OP6], (11) 

where 

, v £ /(<5OS0)d(cOS0) , 

(/> = J ° \i M ' ■ (12) 
Jo d{cos6) 

This normalized function becomes F = when / does not depend on cos#, i.e., in the 
case of a spherically symmetrical distribution. 

First, we present the two test problems. The distributions are calculated for a 
flat space-time and a Schwarzschild space-time. We expect F = 0, for ^ cos 6 ^ 1 
since both space-times have spherical symmetry. We adopt N tot = 2 x 10 9 , r e = 4.5M, 
Too = 500M and v = 0.9 (7 ~ 2) in the numerical calculations. The results are shown 
in Figure HJ There are small irregular deviations from in the number distributions 
F. These small deviations are ~ 0.2% at most, and arise from the finite numbers 
N tot and iVo. In constructing F, we divide the angle into N (= 100) segments. If 
-Wtot(= 2 x 10 9 ) trajectories are randomly distributed, then the mean value is N tot /N , 
and the deviation from the average is 

(iVo/iVtot) 172 ~ 2 x 10 ~ 4 - Additionally, randomly 
chosen emission angles, 9 and 0, are involved in the numerical calculation, so that the 
deviation increases by a certain factor. We numerically checked that the statistical 

— 1/2 

error is roughly proportional to N tot . The error therefore becomes in the limit of 
Ntot 00. Fortunately, this 'noise' level is sufficiently low to allow the effect of black 
hole spin to be examined, as discussed below. 
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Figure 2. The same as 
Figure m but for the extreme 
Kerr spacc-timc. 



(1) • r=4.5M (=re) 

(2) > r=5.5M 

(3) ■ r=7.5M 

(4) ■ r= 10M 

(5) • r=500M 



0.4 0.6 
cos 6 



Figure 3. Evolution of nor- 
malized number distribution 
for different values of the ra- 
dial coordinate r. 



Next, we calculate the number distribution F for the extreme Kerr space-time. 
We adopt the same parameters as those used in the test calculations. The results are 
shown in Figure [2j The number distribution F exhibits a larger deviation, where there 
is a deficit of F ~ —1.7% near the polar regions and an excess of F ~ 0.8% near the 
equatorial region. These values are obviously larger than the statistical error of ~ 0.2% 
in test calculations. The results, however, are opposite to our expectation that the spin 
of the black hole would cause the particles to be collimated along the rotation axis. 

Figure [3] shows how the number distribution F changes with the propagation 
distance. The distributions for some selected radii are calculated from an initially 
uniform distribution. It can easily be seen that the deviation from the average gradually 
increases with r from the emitting radius r e = 4.5M to = 500M. This increase 
almost stops for r > 10M, where the effect of the Kerr space-time becomes very small 
(a/r) 2 ~ 10~ 2 . This result indicates that the particles tend to fall toward the equatorial 
plane due to the Kerr black hole. This feature has been partially discussed by Bicak et 
al.@]. 

In order to express the above geometrical deviation in a simple way, we calculate 
the function F for various values of the Kerr parameter a, the emitting radius r e 
and the emitting velocity 7. It is found that the following empirical relation is well 
fitted to the numerical results of approximately 20 models in the ranges < a < M, 
4.5M <r c < 10M and 2 < 7 < 10 3 , 

F {mB e) ~ -9.2 x 10 (£p (i + ^0 - i) , (13) 

where the angular part is the Legendre polynomial of / = 2. From this relation, the 
most important factor affecting the geometrical deviation is the emission radius r e . Our 
calculation in this section is limited to r c > 4.5M, but the deviation can be increased 
further for particles emitted in the vicinity of the black hole. As shown in Section 
[51 however, the extrapolation of Eq.( ri"3]) for smaller values of r e is not particularly 
meaningful, since the fraction of absorbed trajectories becomes more important. 
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4. Weak gravity with quadrupole moment 

In order to examine the deviation from spherical symmetry, we consider the weak gravity 
limit of Kerr space-time. We introduce asymptotically Cartesian and mass centered to 
order N (ACMC-N) coordinate system^', r', 9', (j)') p]. In the ACMC-2 coordinate 
system, which is the lowest order containing non-spherical terms, the Kerr metric is 
written as follows: 

2M 3Ma 2 cos 2 9' . . 

gw = - 1 + ^ > ( 14 ) 

„, ( 2Masin9' 5Ma 3 sin 9' cos 2 0'\ , . 

gw = r sm9 I — 2 + — 4 j , (15) 



2M AM 2 - a 2 8M 3 - AMa 2 - Ma 2 cos 2 9' 



9W = 1 + — + ~2 + ^3 > ( 16 ) 



= r ' 2 ( 1 + ^ ) 



a 



/2 • 2 nl I 1 , a 



g yy = r' sin 9' 1 



i 2 2Ma 2 sin 2 9' 



.( 2Ma 2 cos9'sm9'\ 
9i>* = r' - 3 J , (19) 

where the relations to the Boyer-Lindquist coordinates (t, r, 9, </>) are 

. g^cos 2 ^' „ „. a 2 cos 9' sin 9' 

r = r'-\ — , 9 = 9 ; , 

2r' ' 2r /2 

= 0', t = t'. (20) 

One of the important non-spherical effects is the quadrupole field in the Newtonian 
level, <7o'0' = — (1 + 2$^), where 

M 3Ma 2 cos 2 9' 

*» = -y+ 2r -3 • (21) 

A particle in this potential field receives an acceleration in the ^'-direction as 
3Ma 2 sin 29' 
~2? 

This is positive in the upper half plane(0 < 9' < 7r/2), while negative for n/2 < 9' < it. 
The quadrupole field therefore bends the orbit of a particle toward the direction of 
equatorial plane. The quadrupole feature of gravity should be one of the mechanisms 
for the anti-collimation shown in Section [31 In the following, we examine numerically 
the magnitude of the deviation in the ensemble of particles which is caused by the 
quadrupole moment. We calculate the particle distribution F(cos#) using Eqs. (jT4l) - 
(TL9l) and geodesic equations, in the same way as done in Section [3] In this space-time, 
we solve the second-order differential equation for the geodesic without using Carter's 
constant Q. The results for a = M,r' e = 4.5M and v = 0.9 (7 ~ 2) are shown in 
Figure HI The number distribution F has a deficit near the polar region and an excess 
near the equatorial region. This behavior is the same as Figure [2j The concentration of 
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Figure 4. Deviation from uniform distribution at infinity for a = M in the ACMC-2 
coordinate, r c = 4.5M and v = 0.9 (7 ~ 2). 
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Figure 5. Escape fraction for various emission radii r e in the extreme Kerr space- 
time. The escape fraction is the ratio of the number of particles which have escaped 
to infinity to the total number of particles emitted from (r c , 6 C ). 



particles toward the equatorial plane can partially be explained by the mass quadrupole 
moment of the central object. The deviation of Figure [2] is roughly twice as large 
as that of Figure HI It is impossible to exactly compare results in different space- 
times exactly. The emission models are different. For example, the constant radius 
r' c = 4.5M in ACMC-2 coordinate does not agree with the radius r e = 4.5M in the 
Boyer-Lindquist coordinate. The correspondence is approximately valid. There is a 
common feature of particle concentration toward the equatorial plane, although exact 
comparison between two different coordinates is impossible, and hence there is still a 
problem in the magnitude. In the weak field limit, the feature is explained by a mass 
quadrupole term. The feature still appears in much strong field, as calculated in Section 

El 



5. Absorption effect 

The number of orbits in which particles are absorbed into the black hole increases as 
the distance between the emitting surface and the black hole decreases. Figure [5] shows 
the escape fraction R cs as a function of cos# e m the extreme Kerr space-time. The 
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Figure 6. Deviation from uniform distribution at infinity. The left panel shows the 
results for emission radius r c = 2.5M, and the right for r c = 1.1M. 

escape fraction is the ratio of the number of escaped particles to the total number 
of particles emitted from (r e ,6 e ). The initial velocity of the particles is chosen as 
v = 0.99999 (7 ~ 224). For a large emission radius, the escape fraction is almost 
unity, which indicates that in-falling orbits to the black hole can be ignored. As the 
radius r e decreases, the escape fraction significantly decreases. Figure [5] shows that the 
distribution of particles escaping to infinity or being absorbed into the black hole is 
not uniform, showing a slight dependence on the angle 9 C . Such dependence indicates 
that the emission distribution can be regarded as being non-uniform. This property is 
important since the non-uniformity at the emission surface naturally leads to the non- 
uniform outflow at infinity. The trajectories of particles emitted at an angle 9 C do not 
necessarily have the same angle 9 = 9 C at infinity, as the correspondence between the 
polar angles is somewhat different at the time of emission and at the time of observation 
due to its modification during the propagation of the particles. It is seldom the case that 
a non-uniform distribution becomes a uniform distribution at infinity, and therefore the 
existence of non-uniformity at emission might be important. 

The results shown in Figure [6] illustrate the effect of the orbits in which particles 
are absorbed in the extreme Kerr black hole. The left and right panels show the results 
obtained for an emission radius r c = 2.5M and r c = 1.1M, respectively. The initial 
velocity is v = 0.99999 (7 ~ 224), and the other parameters are the same as those 
used in the test calculations. Regarding emission for r e = 2.5M, there is an excess of 
F ~ 14% at cos 9 = 1 and a deficit of F ~ —2% at cos 9 = 0. An interesting feature is 
that the excess appears near the polar regions, which is opposite to the result obtained 
for r e > 4.5M, as considered in Section [3j Furthermore, the magnitude is much larger 
than the geometrical effect, which is estimated by extrapolating the empirical formula 
( TT3l) . and is F ~ —6% at the largest deviation point cos 9=1. Thus, it is concluded that 
the number distribution strongly depends on other factors. As shown in Figure El the 
escape fraction for r e = 2.5M varies with a typical amplitude ±10% around the average 
value. This effect overwhelms the modifications which occur during the propagation, 
although the excess at the polar regions is expected to be suppressed to a certain extent 
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if the extrapolation still holds. The right panel of Figure shows the results obtained for 
a smaller emission radius r e = 1.1M. Note that the outer boundary of the ergoregion 
is given by r = M(l + sin#) in the extreme Kerr case, where the emission sphere is 
located within the ergoregion, except for the polar regions. The excess around the polar 
regions transforms into a great deficit of —90%, and the deviation from the average 
also becomes much larger. The origin of this feature can mainly be explained by the 
non-uniformity at the emission surface r e . 

6. Discussion 

A strong gravitational field affects the emission and motion of particles in several ways. 
For example, a central black hole plays an important role in forming the surrounding 
accretion disk, and hence determines the distribution of the particles which are emitted 
into phase space, i.e., the distributions of position and velocity. Strong gravitational 
fields also affect the propagation of particles. In this paper, we have concentrated on 
the latter in a Kerr space-time by calculating the trajectories of particles emitted from 
a 'spherical' shell. Our concern is not the orbit of a single particle, but the behavior 
of an assembly of particles. That is, the deviation from 'spherical symmetry' which 
occurs during the propagation of the particles which are initially uniformly distributed 
on the emission shell and are emitted in isotropic directions. We were able to extract 
the geometrical effect under an ideal situation for the emission at r c > 4.5M. In 
this case, the particles are anti-collimated on average, i.e., they spread in a direction 
perpendicular to the rotation axis in Kerr space-time. This feature can be seen even in 
weak gravitational field with mass quadrupole moment, which partially accounts for the 
physical mechanism. The magnitude depends on the Kerr parameter and the emitting 
radius r e , and the deviation from uniform distribution amounts to a few percent in 
the case of r e = 4.5M. With a further decrease in the radius, the geometrical effect 
increases, but the other effects, such as the absorption of particles in certain orbits, 
become important, and therefore our definition of collimation/anti-collimation becomes 
ambiguous for r e < 4.5M. 

This feature of anti-collimation along the rotation axis does not agree with the 
results of previous works at first sight. The main difference arises from the emission 
models used. Bicak et al. [I] derived orbital parameters of which a particle emitted in the 
radial direction is collimated toward the rotation axis. The condition satisfies Q < 0, 
which corresponds to the 9 motion confined in the range 0<6'<7r/2or7r/2<#<7r 
[TT] . We checked the condition in our numerical simulation, and found that the fraction 
of trajectories satisfying it is typically ~ 0.1%. The isotropic emission is assumed in 
our model, so that most trajectories are irrelevant to the collimation condition derived 
by Bicak et al. [1]. Williamspjl] calculated the production and subsequent propagation 
of e~ e + pairs in the Penrose process within the ergosphere by using a Monte Carlo 
simulation. The results show that escaping trajectories are collimated and that helical 
trajectories exist along the rotation axis. The initial emission within the ergosphere 
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is highly constrained in the Penrose process [TS], since one of the particles in the pair 
should be in a negative energy state. The emission model is therefore different from ours. 
Recently, Bini et al.[16] considered the outflow along the rotation axis, and showed that 
an initially spherical assembly of particles becomes elongated along the rotation axis 
during the propagation of the particles. In their model, the initial direction of the 
velocity is fixed along the rotation axis, which is again different from our model. 

It is not easy to summarize the effects of a strong gravitational field on the 
collimation of emitted particles in a Kerr space-time while taking into account the 
results obtained in other works. In our model, the effect which occurs during the 
propagation has been clarified as a result of simplifying the emission model. We found 
the unanticipated result, i.e., no collimation along the spin axis, but rather spread toward 
equator. This anti-collimation effect along the rotation axis is mainly explained by mass 
quadrupole moment of the Kerr space-time. The geometrical modification integrated 
from the emission point to infinity is ~ (a/M) 2 (6M/r e ) 2,5 , which is not so large. We 
conclude that in order for the collimated outflow to form, the emission model should 
provide a much larger deviation in the emitting direction and in the spatial distribution, 
where the Kerr space-time has a considerable effect. 
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